jordan_sym <-
function(x){
	n <- length(x)
	out <- matrix(rep(0,n^2),n,n)
	out[1,] <- x
	out[,1] <- x
	diag(out) <- x[1]
	return(out)
}

